Resonant transparency of materials with negative permittivity 



o 
o 



E. Fourkal 1 , I. Velchev 1 , C-M. Ma 1 , and A. Smolyakov 2 
1 Department of Radiation Physics, Fox Chase Cancer Center, Philadelphia, PA 19111, U.S.A. and 
Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada 

(Dated: February 2, 2008) 

It is shown that the transparency of opaque material with negative permittivity exhibits resonant 
behavior. The resonance occurs as a result of the excitation of the surface waves at slab bound- 
aries. Dramatic field amplification of the incident evanescent fields at the resonance improves the 
resolution of the the sub- wavelength imaging system (superlens). A finite thickness slab can be 
totally transparent to a p-polarized obliquely incident electromagnetic wave for certain values of the 
incidence angle and wave frequency corresponding to the excitation of the surface modes. At the 
resonance, two evanescent waves have a finite phase shift providing non-zero energy flux through 
the non-transparent region. 
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I. INTRODUCTION 

Propagation of the electromagnetic radiation in ma- 
terials with negative dielectric permittivity (permeabil- 
ity) or the so called left-handed materials (LHM) has at- 
tracted great deal of attention in recent years 0, • 
The increased interest in properties of such media has 
been driven by their potential applications in various 
branches of science and technology. One possible ap- 
plication is related to the possibility of creating the so 
called superlens: a subwavelength optical imaging sys- 
tem without the diffraction limit 0, IE IE 0- Tne su- 
perlens phenomenon is essentially based on amplifica- 
tion of evanescent waves, facilitated by the excitation 
of the surface plasmons 6J. Plasma with overcritical 
density is a simplest example of the negative e mate- 
rial; e = 1 — ojp/uj 2 < for oj < lo v . Phenomena that 
take place in such plasmas are important in a number 
of areas, in particular for the inertial confinement fusion 
(ICF) experiments 0H0. 

In this work we show that the amplification of the 
evanescent waves originating from the interaction of two 
evanescent fields (decaying and growing in space) has a 
resonant character related to the excitation of surface 
modes. Such amplification of the evanescent field al- 
lows the penetration of the electromagnetic radiation to 
depths much greater compared to the incident light wave 
length 0. The evanescent wave incident on the nega- 
tive e slab is strongly amplified when the resonant con- 
ditions for the surface modes are met. For materials 
with the resonance occurs for finite values of 

the wave vector k y (for e = — 1 , the resonance occurs for 
\k y \ — > oo), where k y is the in-plane wave vector com- 
ponent (p-polarization is considered). The presence of 
resonances with finite values of the in-plane wave vector 
k y may significantly improve the overall resolution of the 
subwavelength imaging system. 

Resonant excitation of a surface mode is also the un- 
derlying mechanism behind the absolute transparency 
of a finite thickness slab of material with e < to 
the incident propagating electromagnetic waves (\k y \ < 



w / c ) El ■ In this case, the resonant excitation of sur- 
face modes by the incident light can be achieved via the 
presence of a single transition layer with < e < 1 on one 
side of the slab with e < oEI- The superposition of two 
evanescent waves provides a finite energy flux through the 
region with e < 0, which is equal to that in the incident 
electromagnetic wave. The radiation is then re-emitted 
at the other side of the opaque slab. 



II. INTERFERENCE OF THE EVANESCENT 
WAVES AND EFFECT OF SUPERLENSING 

Non-propagating (evanescent) modes are basic solu- 
tions to the Maxwell's equations for the electromagnetic 
fields in materials with negative dielectric permittivity. 
Often such modes have been neglected assuming that the 
boundary conditions in an infinite medium preclude the 
exponentially growing modes while the decaying modes 
do not contribute to the transmission. It has recently 
been noted however that in a slab of material with neg- 
ative permittivity both (growing and decaying) compo- 
nents are present, resulting in amplification of the evanes- 
cent modes. As shown by Pendry 0, the effect of ampli- 
fication of the evanescent spectrum of the incident light 
(originating from the object) by "negative" materials can 
be advantageously used for subwavelength imaging appli- 
cations, potentially leading to optical system without the 
diffraction limit or superlens. 

The imaging problem can be described in terms of the 
optical transfer function t(x, k yi uj) (k y designates the in- 
plane wave vector directed along the surface of the mate- 
rial), defined as the ratio of Fourier components of image 

field to object field, B^ g (x) / '5^.(0) (for -oo < k y < oo) 
at a given imaging plane x. The transfer function can be 
found by considering a p -polarized wave (electric vec- 
tor in the plane of incidence) incident from vacuum on 
a thin overcritical density slab of thickness d, dielectric 
permittivity e < and magnetic permeability \i = 1 as 
shown in Figure 2] The optical properties of this slab 
are obtained by taking the ratio of the field in the region 
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x > a + d to that at the object plane (in current calcu- 
lations the object plane is assumed to be at x = 0). The 
electromagnetic fields in each region of interest are found 
from solving the well known wave equation, 



dx 



l dB z 
e dx 



0, (1) 



to a coupled surface wave running on opposite sides of the 
slab. As d becomes large so that tanhL/fc 2 , _ ek^d] — > 1, 



(5) 



the two solutions degenerate and we obtain, 
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with a general solution having the following form, 

B z = (A ie lkx + A 2 e~ lkx ) e 1 ^-^ (2) 



where k = uj / 'c^J (e — k 2 c 2 /to 2 ). The electromagnetic 

fields in vacuum regions (x < a and x > a + d) rep- 
resent a sum of incident and reflected waves (x < a), 
and a transmitted wave (x > a + d). Matching solu- 
tions at different boundaries by requiring the continuity 
of B z and l/edB z /dx across interfaces, one arrives at the 
expression for the transfer function, 



T(x,k y ,k ) 



S(d, k y , k ) = 2k 2 e\ 
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A(d, k y , k ) = ((1 + e 2 ) fc 2 - e (1 + e) k 2 ) sin[fcody e - ■£] 

where a; is a distance between the source and the image 
planes and ko = uj/c. Figure [21 shows the absolute value 
of the optical transfer function at a distance x = 2d from 
the source. As one can see, there are in general two sharp 
peaks in the transfer function occurring for certain reso- 
nant values of the in-plane wave vector k y > lo/c where 
the spectrum is completely evanescent. The peaks cor- 
respond to the resonant excitation of surface plasmons 
that are supported by the interface between the thin 
overcritical density medium and vacuum. This means 
that only those fields for which k y > lo/c may resonantly 
couple to this particular surface mode and get dramat- 
ically amplified through constructive interference (these 
are the spectral components that carry sub-diffraction- 
limited resolution of the object). In fact, the zeros of the 
denominator in Eq. [3] 



which is the well known dispersion relation for surface 
plasmons on an overdense medium-vacuum interface. It 
is worth noting here that the transfer function given by 
expression [3] is exact and no approximations were used 
in its derivation. It includes the resonant contribution of 
surface modes with a finite in-plane wave vector k y and 
in the limit e — > — 1, \k y \ >> cj/c it is reduced to that 
obtained in Ref. where the authors investigated the 
resolution limit of the slab of "negative" material with 
dielectric permittivity e = + and magnetic perme- 
ability /i = — 1 + ifi . It should also be noted that the 
location (in k y wave vector space) of the two resonances 
not only depends on the absolute value of the dielectric 
constant e but also the thickness of the slab d. For a 
(3)very thin slab, the two resonances are separated (smaller 
slab thicknesses result in greater peak separation in k y 
space). As the slab thickness increases, the two resonant 
points move toward each other, eventually merging to- 
_ gether forming a single peak. All of these features in the 
k 2 transfer function have not been noticed in the previous 
"rj ] investigations since the authors limited their calculations 
to the asymptotic case k y » lu/c and k y d > 1. As we 
will see below, the presence of such resonances may ad- 
vantageously be used in the design of superlens for im- 
proved image resolution and intensity. 

The transfer function r can be used to find the recon- 
structed field in the image plane in the form 

B img (x,y,t) = I ' S(k v )T(x,k v ,k )e i V'»y- wt ')dky, (6) 



S(d, k y , ko) + A(d, k y , ko) = 

define the dispersion relation for the plasma surface 
eigenmode that is supported by the thin overdense 
medium. The above equation can be reduced to the fol- 
lowing form: 



tanh[W k 2 — ek 2 d] 



2e-i/fco ky 



(l + e 2) fc 2_ e(1 + e)fe 2 



(4) 



In the general case, for a finite value of d there exist two 
solutions to Eq. as seen from Figure @ corresponding 



where S(k y ) is the wave vector spectrum of the source 
(imaged object). Thus, the ability of the system to im- 
age the object is completely determined by the optical 
transfer function which in itself depends on many 
system's physical parameters. In an ideal case, the trans- 
fer function should transfer all spatial harmonics equally 
or t(x, k y , ko) = 1 for — oo < k y < oo. In reality how- 
ever, the transfer function is a non-monotonous function 
of the wave vector k y , medium material type defined by 
e, its thickness d, and the position x of the imaging plane 
relative to the position of the object plane. The super- 
lensing can only be realized for certain set of the above 
parameters. Figure |3] shows the optical transfer function 
for two layers of materials with dielectric permittivities 
e x = -l+0.001iand e 2 = -1.0292+0.001i at the imaging 
plane x — 2d. As one can see the transmission band (in 
the wave vector space) is significantly broadened in the 
second material (e = —1.0292 + O.OOli) due to the reso- 
nant excitation of a surface plasmon with a finite value of 
the wave vector fej«4 Such surface plasma wave is not 
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excited in material with e = — 1 (assuming the slab thick- 
ness d such that the condition tanh[^/ — ek 2 d] = 1 is 

satisfied) , as can be seen from the dispersion relation |JSJ • 
Let us estimate how well we can image an object using 
these two materials. For a sake of simplicity we assume 
that our object is represented by two slits of a certain 
width located at a given distance away from each other 
as shown in Figure ^ Substituting the Fourier trans- 
form of the object together with the transfer function (J3J) 
into Eq.0we arrive at the reconstructed image shown in 
Figure for the case when the incident light wavelength 
is A = 350 nm. As one can see, both materials pro- 
vided considerable focusing, yielding the sub- wavelength 
resolution of the object. However, as we expected the 
image resolution and its intensity for the second layer 
(e = —1.0292) is superior to that using material with 
e = —1. This suggests that there is a range of sys- 
tem parameters for which a significant focusing can be 
achieved 4] necessitating a further parametric study in 
order to understand the relation between different phys- 
ical parameters of the system. 



III. SURFACE WAVE INDUCED TOTAL 
TRANSPARENCY OF MATERIAL WITH 
NEGATIVE PERMITTIVITY 

Surface wave induced amplification of incident elec- 
tromagnetic waves with spectral components of in-plane 
wave vectors k y satisfying the condition \k y \ > oj/c has 
been considered in the previous section. It was shown 
that the resonant amplification occurs as a result of the 
excitation of a surface wave. A slab of negative e ma- 
terial surrounded by vacuum supports such a surface 
mode for which its phase velocity is always sub-luminal 
or u/ky < c. As a result, only those modes of the inci- 
dent light for which \k y \ > uj/c have been amplified. It is 
possible however to amplify the propagating modes with 
\k y \ < oj/c too, thus creating the conditions for the abso- 
lute transparency to the incident propagating wave. This 
can be done by creating conditions for the excitation of 
a surface mode with phase velocity greater than that of 
the speed of light. 

Consider p-polarized light obliquely incident from vac- 
uum on a two-layer structure having dielectric permittiv- 
ity distribution shown in Fig. Such a system can be 
formed by placing an undercritical density plasma layer 
(with thickness d and electron density corresponding to 
the plasma frequency oj p ,i) to an overcritical density 
plasma layer (with thickness a and electron density cor- 
responding to the plasma frequency w P ,2)- The plasma- 
plasma interface supports a surface wave with disper- 
sion relation found from the solution to the Maxwell's 
equations [J , 



where a 2 = fc 2 - e^ 2 /c 2 , e; = 1 - uJ 2 :i /to 2 , (i = 1,2); 
A = w 2 i/Wp ;2 ; w is normalized to the plasma frequency 
in the region where e < and k y is normalized to the 
classical skin depth S = c/lu Pi 2- It can be easily seen 
that the phase velocity of the surface wave on a plasma- 
plasma interface can be greater than that of light, so that 
they can couple to radiating electromagnetic fields. This 
means that the incident p-polarized electromagnetic wave 
may excite a surface mode on a plasma-plasma interface 
if the resonant condition (external field frequency and its 
wave vector's tangential component have to match those 
that are determined from the dispersion relation for the 
plasma surface wave) is satisfied. The optical properties 
of this dual-layer system can be found from matching 
the fields of the incident/reflected electromagnetic waves 
with those of the surface wave at the three interfaces. 
The electromagnetic field in each region B = (0,0, B z ) 
is a solution to the wave equation ^ with general solu- 
tion given by expression^ The electromagnetic fields in 
vacuum regions (x < — d and x > a) represent a sum of 
incident and reflected wave (x < —d), and a transmitted 
wave (x > a). Matching solutions at different boundaries 
by requiring continuity of B z and 1 /edB z /dx across inter- 
faces, we obtain the unknown expansion coefficients with 
the transmission coefficient having the following form, 



T 



4e i(fc 1 d+* :3 a-(c t +(i)fcoco8e) fcofclfc2eie2 cos g 



(8) 



ai/ei + a 2 /e2 = 



VI + A - 2uj 2 ' 



(J) 



(H-T) 

3? = (1 + e 2lkld ) fc fcifeeie 2 cos 6 + e 2lk2Cl (l + e 2lkld ) 
kQk\ki(.\t2 cos# 

T = 2ie^ kld+k2a) k 2 e 2 cos[fc 2 a] (kf + fc 2 e 2 cos 2 6) sin[fcid] 
+feiei cos[fcid] {k 2 + fco e 2 cos2 sin[fc2a] 
— iko (fc 2 e 2 + k 2 e 2 ) cos 8 sin[fcid] sinffeaa]. 

where 9 = axcsin[k v /ko\. Figure0shows the transmission 
coefficient as a function of the incidence angle. As one 
can see there is a sharp increase in the transmission prop- 
erties of the system when the angle of incidence matches 
certain resonant value of 6 at which point the transmis- 
sion coefficient reaches unity. The resonant value exactly 
corresponds to that given by an expression J7J for the dis- 
persion relation of the surface plasma wave. Thus, the 
anomalous transmission occurs even for a system consist- 
ing of an undercritical density plasma layer adjacent to an 
overcritical density plasma slab, so that there is no need 
to form a sandwich-like structure as argued in Refs. 

0113 

in which a second surface plasma wave has to be excited 
on the opposite side of the overcritical density layer to 
achieve the same effect. In other words, the anomalous 
light transmission can be achieved through excitation of a 
single surface plasma wave. The transmission coefficient 
is also a non-monotonous function of both plasma layer 
thickness with a single maximum reached at certain cor- 
related values determined from the interference condition 
between evanescent fields in both plasma slabs. 
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IV. EVANESCENT WAVE INTERFERENCE 
AND THE ENERGY TRANSPORT 

It is often assumed that the evanescent waves do not 
carry the energy. Therefore, in the problem of total trans- 
parency of a layered structure there occurs a question of 
how the energy is carried through the non-transparent 
media where the only solutions are the evanescent modes. 
One must remember however that the general solution in- 
side the negative e medium is a sum of two exponential 
functions, one that decays with the distance ~ e~ x [x 
points in the propagation direction) and the other grows 
~ e x . For the superposition of decaying and growing 
modes, E,B ~ A\ exp[— ikx] + A 2 exp[ikx] (k is purely 
imaginary decay constant) , the x component of the time 
averaged Poynting vector S x 

S x = ^Re[E y B*] ~ Re[k(A 1 A 2 - A 2 A\)] ~ 2Im[A x AZ\, 

may become finite when the combination A\A* 2 has a fi- 
nite imaginary part, which requires a finite phase shift 
between A\ and A 2 . Therefore, a finite energy flux oc- 
curs as a result of the superposition of two evanescent 
modes with a finite phase shift . We shall call this the 
interference of the evanescent modes. 

It is easy to show that the required phase shift can 
be obtained when two evanescent modes inside the neg- 
ative e region are matched to the outgoing (transmitted) 



wave in vacuum. Matching of the evanescent solutions 
with the incident vacuum wave (at the other side of the 
negative e region) shows that the total transmission may 
be obtained only when the transition layer with e > is 
included|ll|. The condition of the absolute transparency 
is equivalent to the resonant condition for the excita- 
tion of the surface plasma mode. At the resonance the 
Poynting flux inside the slab becomes equal to that of the 
incident radiation and the opaque plasma slab becomes 
absolutely transparent. This can be a possible mecha- 
nism of the anomalously high transparency of overdense 
plasma recently observed in the experiment [T(| . 



V. SUMMARY 

In conclusion, we have shown that the excitation of 
surface modes leads to the resonant transparency of opti- 
cally opaque materials. Presence of such resonances may 
improve the resolution and the signal intensity of the sub 
wavelength imaging system. It has been also shown that 
the resonant excitation of surface modes is an underly- 
ing mechanism behind the total transparency of an over- 
dence plasma slab to the incident electromagnetic wave. 
The energy flux through the negative e region occurs as 
a result of the interference of the evanescent modes. 

This work is in part supported by Strawbridge Family 
Foundation, and NSERC Canada. 
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FIG 1: Schematic geometry of the dielectric constant distri- 
bution for the case of a single planar medium. 
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FIG. 5: Field distribution in the image plane for two slabs 
with ei = -1.0292 + O.OOli (solid line) and e 2 = -1 + O.OOli 
(dotted line) and thickness d = 1.8c/u p . Magnetic permeabil- 
ity of both materials is assumed to be {i = 1. The wavelength 
of the incident light is A == 350 nm. 



FIG. 2: The absolute value of the optical transfer function 
at a distance x = 2d from the source as a function of the in- 
plane wave vector k y . The dielectric constant of the medium 
and its thickness are e = —1.0292 + O.OOOOli and d — 1.8c/u p 
correspondingly. k y is normalized to classical skin depth c/u p . 
External radiation frequency ko = 0.702. 
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FIG. 3: The absolute value of the optical transfer functic 
at a distance x = 2d from the source as a function of the h 
plane wave vector k y for two different materials with dielectr 
permittivities ei = —1.0292 + O.OOli (solid line) and £2 
— 1 + O.OOli (dashed line) and thickness d = 1.8c/ Up. k y 
normalized to classical skin depth c/lo v 



FIG. 6: Schematic diagram of the spatial density (dielectric 
constant) distribution for the case of a layered slab. 
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FIG. 4: The field distribution in the object plane. The ob- 
ject comprises two slabs of thickness 15 nm, separated by a 
distance 80 nm from their centers. 



FIG. 7: Transmission coefficient as a function of the inci- 
dence angle. A double-layer plasma system becomes com- 
pletely transparent at the incidence angle #=0.671955. Un- 
dercritical and overcritical plasma slabs dielectric constants 
and thicknesses are ei = 0.3428, d = 27 * c/u P 2 and £2=- 
2.97, a = 3.12 * c/ui P 2 correspondingly. External radiation 
frequency u/u P 2 = 0.5019. 



